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Abstract 

Bivariate generating functions for various subsets of the class of per- 
mutations containing no descending sequence of length three or more are 
determined. The notion of absolute indecomposability of a permutation 
is introduced, and used in enumerating permutations which have a block 
structure avoiding 321, and whose blocks also have such structure (recur- 
sively). Generalizations of these results are discussed. 



1 Introduction 

The association of a permutation of {1, 2, . . . , n} with its graph provides a ge- 
ometric viewpoint in which to consider pattern avoidance. Thus, for example 
the permutation avoids 132 if the points to the left of the highest point, all lie 
above the points to the right, and this condition is true recursively of the points 
to the left and to the right of the highest point. This condition is illustrated in 
Figure |l|. Likewise any set of n points in the plane, no two of which lie on a 
horizontal or vertical line can be associated with a permutation of {1, 2, . . . , n}. 
In that case we can represent the geometric information about 132 avoidance 
as a picture, with a point for the maximum, and two bounding rectangles to 
its left and right, the former lying above the latter, together with an implicit 
understanding that the structure within the rectangles is to be similar. These 
associations provide a geometrical context in which to consider pattern avoid- 
ance, and are a common tool in understanding the class of permutations that 
avoid one or more patterns. We explore some ramifications of this viewpoint 
in the very simple situation of 321 avoiding permutations (which we denote by 
A(321)). 

The main purpose of this paper then, is to show how the geometric context 
provides a simple method to obtain more detailed enumeration results about 
321 avoiding permutations than have hitherto been available. Moreover, these 
results are obtained uniformly in some sense. The underlying technique con- 
sists of identifying a suitable geometric configuration which must be attained or 
avoided, and then using the structural constraints which that implies in order 
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Figure 1: Schematic representation of a 132 avoiding permutation. 



to compute the generating function, often multivariate, of the associated collec- 
tion of permutations. This technique has applications beyond the scope of 321 
avoidance, a few of which we consider in the final sections of the paper. 

Thought of as a set of points, a permutation avoids 321 if it does not contain 
three points, every pair of which determines a line segment of negative slope. Of 
course it is also the case that any 321 avoiding permutation is the merge of two 
increasing sequences. It is easy to see that one of these sequences can be taken 
as the sequence of left to right maxima, that is, those elements which dominate 
all of their predecessors. These elements now determine a sequence of vertical 
and horizontal ranges in which the remaining elements of the permutation must 
lie. The situation is illustrated in Figure [l]. Subject to having a fixed number of 
left to right maxima, the possible 321 avoiding permutations that can be formed 
are in one to one correspondence with assignments of non-negative integers to 
the cells of that diagram, so that no two cells whose centres are connected 
by a segment of negative slope are assigned positive labels. We refer to this 
diagram, with the cells containing positive labels simply marked (but the labels 
themselves suppressed) as the skeleton of the permutation. 

As is well known, the total number of 132 avoiding permutations of length n and 
the total number of 321 avoiding permutations of length n are the same, both 
being equal to the nth Catalan number. However, the schematic representation 
of the 132 avoiding permutations makes the correspondence between them and 
plane binary trees clear and hence also the equation satisfied by the generating 
function of the class, while the corresponding diagram for the 321 avoiding 
permutations does not. In some sense ^4(321) is a class which exhibits more 
subtle structure than ^4(132) does. 

Alternatively we might replace each cell in Figure [ljwith its central point, pro- 
ducing a triangular subset of the integer grid. Given a 321 avoiding permutation, 
the points corresponding to cells with positive labels determine a path in this 
grid, lying on or below a diagonal line and having steps which are of the form 
(a, b) where a and b are non-negative integers, not both 0. This path is also 
shown in Figure [ij Such paths are enumerated in jl5| and H], the former of 
which also provides a bijective interpretation of the relationship between their 
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Figure 2: A 321-avoiding permutation with four left to right maxima. The given 
occupancies of the cells define the permutation 45127361089. Note that the 
relative horizontal position of the occupants of the two vertically aligned cells 
is determined by the fact that the resulting permutation is not to contain 321. 

enumeration and that of standard Dclannoy and Schroder paths. 

Our objective is to find the enumeration via generating functions of various 
subsets of the 321-avoiding permutations. Generally we will aim to obtain a 
multivariate generating function (in fact, at most bivariate) say in x and y 
where the coefficient of x n y k might be the number of 321-avoiding permutations 
of a certain type having n left to right maxima and k other elements. Then 
simple substitutions will allow us to compute either the associated univariate 
generating functions, or permutations of size n having k left to right maxima, or 
other similar variations. Whenever possible, the name of a generating function 
will reflect the name of the class that it enumerates, so that for example the 
univariate generating function for the class A(321) will be A(i), while a bivariate 
form might be A(x,y). 

To carry out the enumerations of various structurally defined subsets of the 
collection of all 321 avoiding permutations we will make use of the structural 
relationships that hold between the subsets and the class as a whole, and the 
corresponding algebraic relationships which hold for the generating functions. 
That is, our methodology is firmly in the school represented by Jl0|], or ||. 

In the next section we rederive equations describing the basic generating func- 
tions for 321 avoiding permutations by defining a context-free language whose 
elements are in one to one correspondence with 321 avoiding permutations. This 
method is used because it is no longer than extending the results of jig] to the 
bivariate case, and also because it represents a technique of wider applicability 
in the field of pattern class enumeration. 

The following section then applies the results to the problem of enumerating the 
subsets of 321 avoiding permutations consisting of: plus irreducible, minus irre- 
ducible, plus indecomposable, and absolutely irreducible permutations. These 
terms, and their significance for enumeration questions are defined below. In 
the univariate case, only the last of these is definitely new, although we have 
not found detailed expositions of the others in the literature. 
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Figure 3: An illustration of the correspondence between occupied cells in the 
triangular grid and lattice paths. Points in the subset are marked as solid circles, 
and the corresponding path is illustrated. The associated word is mhumumhdd. 

In the penultimate section we make use of the enumeration of the absolutely 
irreducible 321 avoiding permutations to enumerate another class, built from 
321 avoiding permutations by a recursive construction based upon the wreath 
product introduced in ||. 

In the final section we try to foreshadow future applications of the methods 
illustrated here, and mention some connections with other work in the area of 
pattern classes. 

2 Enumeration of marked paths 

In Figure || we rotate the triangular grid from Figure [l] clockwise by 45°, and 
then reflect it through the x-axis. The paths we then obtain begin from (0,0), 
and travel through lattice points (the sum of whose coordinates is even) along 
segments of slope lying between —1 and 1 inclusive. We will introduce a gram- 
mar which describes such paths. The purpose of transforming the original dia- 
gram is purely psychological. The grammar which we describe is the language 
accepted by a certain deterministic pushdown automaton, and the lattice path 
then models the number of elements held in the stack of this automaton as the 
word is processed. Given our overlying motivation of considering pattern classes 
this also provides a nice link to the generation of 312-avoiding permutations us- 
ing a stack. 

Allowing path segments of arbitrary rational slope in [—1, 1] would obviously 
be a problem in a finite language. We avoid this problem by adding marks to 
the end of each path segment, and then replacing each segment by a horizontal 
segment followed by one at an inclination of ±45°. This is illustrated in Figure 

I 

We now describe a context free grammar which generates a language that de- 
scribes all, and only, the lattice paths which correspond to marked subsets of 
the triangular grid of size n for some n. Constructing this grammar is relatively 
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straightforward. We use four terminal symbols, three of which stand for unit 
segments in the path, and one of which represents a vertex in the set: h for 
horizontal segments, u for upward segments, d for downward segments, and m 
for marking a vertex of the original path. 

Informally words in the language are described as follows: 

• Any change of direction after a d or a u requires an intervening m. 

• h's may only occur in blocks immediately following an m or at the begin- 
ning of a word. 

• In any initial segment there must be at least as many u's as d's, but in 
the whole word the total number of each is the same. 

• m's can occur anywhere. 

This description is easily formatted into a grammar, modelled on the standard 
grammar corresponding to Dyck paths (which do not allow horizontal steps, and 
do not require marks). Each non-terminal symbol in the grammar represents an 
excursion, that is a path starting and ending at the same level and not passing 
below that level. These excursions are sometimes restricted by the immediately 
preceding symbol. 

S -> e\hS\ uUdD \ mM 

U -> mM | uUdD 

D -> e | mM 

M -> e\\\S\ uUAD 

As in H each non-terminal symbol of the grammar is associated with a gener- 
ating function (denoted by the same symbol) in variables h,u,d,m. This gen- 
erating function is obtained by taking the sum of the monomials corresponding 
to words represented by that non-terminal. 

The grammar above is clearly unambiguous since in each rule the initial symbols 
of differing productions differ from one another. So, it is a simple matter to 
obtain a system of equations satisfied by the generating functions of the non- 
terminals, namely: 

S = l + hS + udUD + mM 
U = mM + udUD 
D = 1 + mM 
M = l + hS + udUD. 

Using a symbolic algebra package or, in a pinch, by hand, this system can be 
solved. We are interested principally in the function S describing words of the 
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language, and this is described as S = (1 + m)Si where: 

(—udm — udm + udmh + udm h) Sf + 
(1 — h — mh — ud + udh + udmh) Si — 1 + ud =0. 

Given a word in the language, the number of left-to-right maxima in the 321 
avoiding permutation of which it is the skeleton is equal to one more than the 
sum of the number of h's and the number of u's. Also, the size of the set that 
it encodes is equal to the number of m's. So, we can reduce to the generat- 
ing function S(x,y) where the coefficient of x n y k is the number of skeletons 
with n left to right maxima, and k internal marked cells through the following 
substitutions: 

h — ► x, u — > x, d — > 1, m — * y, 

followed by multiplication by x (and addition of 1 for the empty graph). This 
yields: 

S(x,y) = l + x(l + y)S 1 (x,y) 

= (xy + xy 2 )Sf + (xy + x - 1)S\ + 1. 

The first of these equation can be solved for S\ , with the result being substituted 
in the second. After some further simplification this yields: 

xS 2 + {xy + x-2y-l)S+l + y = 0. (1) 

By substituting y = 1 we will obtain the total number of allowed markings. So, 
defining S(x) = S(x, 1): 

S^x) 2 + (2x - 3)S(x) + 2 = 0. 

The discriminant of the latter equation is v '4x 2 — 12.x + 1 illustrating a connec- 
tion between these numbers and the Schroder numbers (sequence A001003 of 
fU). In fact: 

S(x) = 1 + 2x + J2 2 n+1 s n x n 

n=l 

where s n is the nth Schroder number. 

This sequence of coefficients also arises as the number of non- crossing graphs, 
that is, graphs with n vertices arranged as the vertices of a convex polygon, 
with straight edges connecting these vertices subject to the condition that no 
two edges should intersect at an interior point. This result is due to 0, with a 
more modern derivation, as well as other related results, given in ||. 

If we consider the ("J 1 ) cells of the original grid as the vertices of a graph, T n , 
two vertices being adjacent if they are connected by a line of negative slope, 
then the coefficient of x n in S(x) counts the independent subsets of T n . The 
number of non-crossing graphs is also the number of independent sets in a graph. 
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Namely, take as vertices of the graph the possible edges two such vertices being 
adjacent if the segments which they represent meet internally. In this graph 
NC n , a non-crossing graph corresponds to an independent set. 

So, the generating function for independent subsets of the sequence of graphs T n 
and NC n +i are the same. In fact, inspection of the results in |)j together with 
a little algebra shows that this is also true of the bivariate generating functions 
which mark the sizes of the independent subsets. That is: 

Proposition 1 For every n and every k, T n and NC n +i have exactly the same 
number of independent subsets of size k. 

However, it is easy to see that for n > 4, T n and NC n +\ are not isomorphic. For, 
T n has exactly four isolated vertices, while NC n +\ has n + 1 isolated vertices. 
Detailed expressions for the coefficients in S(x) and S(x, y) can be found in Q 
(Theorem 2, part (ii)) as well as discussions of their asymptotic expansions. 



3 Consequences for 321 avoiding permutations 

Before turning to the enumeration of various subsets of the 321 avoiding per- 
mutations we begin with some remarks about the full class (whose enumeration 
is, of course, already well understood beginning apparently from |TT|). 

In our original setting, the marked cells arose by considering a 321 avoiding per- 
mutation having n left-to-right maxima. We argued that any such permutation 
corresponded to a labelling of marked cells with positive integers representing 
the number of elements of a permutation contained in a particular cell. Let 
A(x,y) be the generating function for 321 avoiding permutations where the 
exponent of x denotes the number of left to right maxima, and that of y the 
number of remaining elements. Since we obtain a 321 avoiding permutation from 
its skeleton by replacing a single cell, marked by a y in S(x,y), by a positive 
integer, marked therefore by y n for some n > 0, we obtain: 

A{x,y) = S(x,y + y 2 H ) = S[x, ' 
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We can also make this substitution in the equation that S satisfies and then 
simplify to obtain: 

yA 2 + (x - y - I) A + 1 = 0. (2) 

On the other hand, it is perhaps more natural to count permutations of a com- 
mon size. So, using A am to denote the generating function where the coefficient 
of x n y k is the number of 321 avoiding permutations of length n having k left to 
right maxima, we obtain: 

A am (x,y) = S (xy, — — 
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By algebraic manipulation this function also satisfies a quadratic equation with 
coefficients polynomial in x and y namely: 

xA 2 am + (xy-x- l)A am + 1 = (3) 

A further reduction in complexity occurs when we substitute y — 1 in A am (or 
V = x in (§)) giving: 

xA{x) 2 - A(x) + 1 = 

thus confirming, in a rather roundabout way, that the total number of 321 
avoiding permutations of length n is enumerated by the Catalan numbers. 

The coefficient of x n y k in A am (x, y), which is non-zero only for 1 < k < n is a 
Narayana number. 

These numbers also arise in but not as a direct translation of this result 
since we are no longer in the context of path counting. They also arise in a 
number of other contexts including the enumeration of A:- way trees ((|]) and as 
the number of non-crossing partitions of n (0). 

We now turn to the enumeration of various subsets of A(32l). First let us define 
those classes and the symbols used to specify them: 

Definition 2 Let n be a permutation (in A(321)). Then: 

(A + i rr ) 7r is plus irreducible if it does not contain a subword of the form i(i+l), 

(A-i rr ) 7r is minus irreducible if it does not contain a subword of the form 
*(*-!), 

(A + i nc i) 7r is plus indecomposable if it does not have a proper initial segment 
whose values form an initial segment of [1, n], 

(A-i n d) 7r is minus indecomposable if it does not have a proper final segment 
whose values form an initial segment of [1, n], 

(Ai rr ) it is absolutely irreducible if it does not have a proper subword of length 
greater than 1 whose values form an interval in [1, n]. 

The irreducible or indecomposable elements of a collection of permutations can 
(under suitable closure properties) be thought of as components in the con- 
struction of the other elements of that class. Again, granted certain closure 
and uniqueness assumptions, this can allow enumeration of the entire set based 
on an enumeration of one of the collections of components, or vice versa. This 
particular exposition of a general combinatorial theme is explored in We 
note that the results in that paper could be used to derive the univariate gen- 
erating function for the plus irreducibles and plus indecomposables in A(321) 
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(results which we will rederive here as a result of obtaining the bivariate form) . 
Furthermore, the only minus decomposable permutations that avoid 321 are of 
the form: 

(fc + l)(/c + 2)---nl2---fc 
so we will not concern ourselves with that case. 

The condition of absolute irreducibility is a new one, and we will see its appli- 
cation in the next section. The definition is not so unnatural as it might appear 
to be at first sight. In terms of the graph of a permutation it says that if some 
proper, non-singleton, part of the permutation is bounded by a rectangle, then 
there must be at least one element of the permutation outside of the rectangle 
but in either the vertical strip or the horizontal strip determined by it. 

Enumeration results in this section generally take equation ([!]) as their starting 
point. Recall that this provides the generating function S(x,y) for skeletons 
of 321 avoiding permutations, with the exponent of x marking the number of 
left to right maximals, and that of y the number of occupied cells. So, all the 
generating functions we compute will be in the form where the coefficient of 
x n y k marks the number of permutations of that type having n left to right 
maxima and k other elements. As usual, a simple change of variable, replacing 
x by xy and y by x would produce the function enumerating by total number 
of elements, and number of left to right maxima. 

If 7r is a plus irreducible member of A(321) then no cell can be occupied by 
more than one element. Among the diagrams that meet this criteria, the plus 
reducible elements contain sequences of more than one left to right maximum 
such that the vertical and horizontal bands which they determine are otherwise 
empty Suppose then that we knew the generating function A + i rr (x,y) for the 
plus irreducible members of the class. The preceding sentences imply that we 
would obtain the generating function S(x,y) by replacing x in A +irr (x,y) by 
x/(l — x). So, since the inverse of sending x to x/(l— x) is to send it to x/(l+x): 

A +irr (x, y) = S(x/ (1 + x),y). 

Substitution and simplification in equation |l| then yields: 

x(y + l)A\ lrr - {xy + 2y + l)A +lrr + (x + l)(y + 1) = 0. (4) 

The corresponding univariate form is: 

x[x + l)A 2 +lrr -(x+ l) 2 A +vrr + (x + l) 2 = 0. 

An element 7r of A(321) can only be minus reducible if some left to right maxi- 
mum k is followed immediately by k — 1. So 

7r = ak (k — 1) /3, 

for some a, (3 £ A(321) (with j3 of course having all its values increased by k). 
We can make this decomposition unique by requiring k, k— 1 to be the first pair 
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of elements witnessing minus reducibility. Then a is minus irreducible, while [3 
could be any 321-avoiding permutation. Thus we obtain: 

A(x,y)=A 

—irr 

(x, y) + A 

—irr 

(x,y)(xy)A(x,y). 

Or, solving for A^i rr (x, y): 

A— i rr (x, y ) = 

The bivariate algebraic equation for A 
form is more presentable: 

(X 4 + X 2 + x)A 2 _ lrr + (1 - 2x 2 )A^rr + 1=0. 

Enumerating plus indecomposables is easier and standard. Every element of 
^4(321) is either of length or of the form a\ai ■ ■ • ot c where each a% is a plus 
indecomposable, shifted upwards by the sum of the lengths of the preceding 
a's. Since this decomposition is unique, then using A +in d to enumerate the non 
empty plus indecomposables, we obtain: 



1 -^-+ind 

which can then be readily solved for A + i nd . 

Finally we come to absolute irreducibility. Since the absolutely irreduciblcs 
form a subset of the collection of plus indecomposables, and of the plus irre- 
ducibles, we begin with the form of the skeleton function which is like that for 
plus indecomposables. This already reduces us to permutations that are plus 
indecomposable, and plus irreducible in their non left to right maxima. 

b +ind (x,y) = r. 

1 + S(x,y) 

Which, by now standard manipulations, satisfies: 

(1 + v)S 2 +md - (1 + x + xy)S +md + xy + x = 0. 

Consider which non empty rectangles in the diagram associated to an element of 
S+i n d might not contain other elements inside the vertical and horizontal strip 
which they define. In order for this to hold, the top edge of the rectangle cannot 
cross a vertical line in the triangular grid of cells, nor can the left edge cross 
such a horizontal line. So, the upper right and lower left corners lie outside of 
the grid. Such a rectangle is illustrated in Figure |[ The vertical area above 
the rectangle is automatically empty as is the horizontal area to the left. So 
problems can occur only when we have a non-empty sequence of left to right 
maxima such that there are no marked cells in the horizontal or vertical strip 
which they define. 



A ( x --y) (5) 

l + xyA(x,y)' 

.irr is not very pretty, but the univariate 
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Figure 4: A potential rectangle for the violation of indecomposability. 



If we knew the function Ai rr (x, y) how could we compute S+i n d(x,y)7 An 
element of the latter class could be obtained beginning from an element of 
Ai rr (x,y) by inflating some of the left to right maxima into a sequence of such 
maxima, adding no additional elements in the horizontal or vertical strips which 
they determine. If we imagine in Figure ^ that the illustrated rectangle (and 
subrectangles of it) are the only ones which cause a violation of absolute in- 
decomposability, then that permutation has been constructed by inflating the 
left to right maximum just to the left of the rectangle into six such maxima. 
As we've already insisted on plus indecomposability, a rectangle whose leftmost 
boundary is to the left of the first maximal cannot be problematic, and so there 
always is an available maximal to inflate. 

There is just a single exception. The permutation 1 is absolutely indecompos- 
able, but when we inflate it we do not obtain plus indecomposable permutations. 

Thus, beginning from Ai rr (x, y) — x, we should replace x by xj (1 — x) in order 
to obtain S + i n d(x,y). Inverting this replacement we get: 

A irr {x,y)=S +ind (^- x ^-^. 

Carrying out these substitutions and manipulations on the equations satisfied 
by the generating function yields: 

(1 + x){l + y)A 2 rr + (xy - l)A irr + xy = Q. (6) 

Or in univariate form: 

(a; + l) 2 A 2 rr + (x 2 ~ l)A irr + x 2 = 0. 

In fact, this does not quite get us all the irreducibles as it omits 1, 12, and 
the empty permutation. Adjusting the equations to include this one comes at 
considerable cost to their appearance, so we prefer to leave the equation as it 
stands, adding the necessary 1 + x + x 2 to the generating function post facto. 
Table || summarizes the sizes of these subsets of ^4(321). 
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n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


all 


1 


1 


2 


5 


14 


42 


132 


429 


1430 


4862 


16796 


plus irr. 


1 


1 


1 


2 


4 


9 


21 


51 


127 


323 


835 


minus irr. 


1 


1 


1 


3 


10 


31 


98 


321 


1078 


3686 


12789 


abs. irr. 


1 


1 


2 





2 


2 


7 


14 


37 


90 


233 



Table 1: Sizes of A(321) and some of its subsets 



4 A "fractal" class 

As an application of the final results of the preceding section we will show how 
the knowledge of the generating function for A irr can be used to compute that 
of a much more complicated class. At the risk of further abusing a term which 
has suffered much abuse already we would like to introduce the class F(321) 
of fractal 321-avoiders. These are permutations which, from a distance, appear 
to avoid 321 but which on closer inspection are made up of blocks, arranged in 
a 321-avoiding pattern where each block appears to avoid 321 but perhaps on 
closer inspection is in fact made up out of blocks . . . 

That is, 7r e ^(321) if cither, n e A(321), or n = ci\U2 ■ ■ ■ ot c where 

• the values occurring in each on form an interval, 

• the permutation represented by on is in F(321), and 

• the relative ordering of the on , interpreted as a permutation of length c is 
in A(321). 

Many well-known permutation classes can be defined as fractal classes in this 
way, or occasionally as natural subclasses of such fractal classes. For example, 
the class of separable permutations is precisely the fractal class generated from 
the finite base class {12, 21}. 

There is a complementary bottom up description of F(321). Namely, this class 
is the closure of the class consisting just of 1 under the operation of replacing 
an element of a permutation by a 321-avoiding block. For example: 

1 -> 2413 -> 423615 

where we initially replace 1 by 2413 and then replace the element 2 by the 
permutation 312 while retaining its relative order within the entire permutation. 
Geometrically, we begin with the graph of 2413 and then expand the vertex 
representing 2 into a copy of the graph of 312. Such replacements could just 
as easily be applied to each element of a permutation and, in some sense, they 
already have been, only 1 has been replaced by 1 in three instances. Thus 
the two descriptions arc equivalent the permutation 423615 consists of blocks 
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(423) (6) (1) (5) whose relative order is 321-avoiding, and where each block is in 
F(321) (in this instance, in fact in .4(321)). 

We could also define -F(321) algebraically using the wreath product operator of 
[|| as the smallest non-empty class X satisfying the equation X = X I A where 
A = 4(321). This corresponds to the bottom up description, while the top 
down one would suggest X = A\X. Consider the first equational description of 
X. Since X contains A we also get that X contains A\A, Then also X contains 
(A l A) l A and so on. Letting A n = A n ~ x lAfbrn>l, and A 1 — A we obtain 

X D u~ A». 

On the other hand, the right hand side is contained in its wreath product with 
A, and so by the definition of X: 

X = U~ A n . 

The second equational definition can be manipulated in the same way and in 
fact, as the wreath product is associative, leads to the same equation, thus 
confirming that the two approaches are indeed equivalent. 

Such an algebraic representation suggests that we ought to be able to transfer 
our knowledge of generating functions for A to similar knowledge about F. 
There is though, a small complication. This arises from the fact that the choice 
of blocks to witness the fact that a permutation belongs to F(321) is not uniquely 
defined. We need to obtain uniqueness of some sort if we hope to carry out the 
enumeration, and the following general result helps to provide that. 

Definition 3 Let 9 be a permutation of length k and let ir be a permutation 
of length n. Then tt is ^-decomposable if n = a±a2 • ■ • ot-k f 0T some non-empty 
subwords oti such that the set of values occurring in each of the on forms an 
interval and the relative ordering of these values agrees with the relative ordering 
of the corresponding elements of 9. The factorization tt = ct\ct2 ■ • • O-h is called 
a ^-decomposition of tt. 

With this new definition, we see that a permutation of length 3 or more is plus 
decomposable if and only if it is 12-decomposable, while a permutation tt is 
absolutely irreducible if and only if it is not ^-decomposable for any 9 ^ ir. 

Proposition 4 Let n be an arbitrary permutation. Then there is a unigue 
absolutely irreducible permutation 9 such that tt is 9 -decomposable. Moreover, 
if 9 =/= 12 and 9 ^ 21 then the 9 -decomposition of tt is also unique. 

For example, for 423615 this decomposition is (423)(6)(1)(5) with relative or- 
dering 2413, while for 724513986 it is (7) (24513) (98) (6) with relative ordering 
3142. On the other hand 123 which is 12-decomposable admits two such de- 
compositions. 
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Proof: Let tt be given, say tt = p\P2 • • -Pn- To each pi associate a maximal 
proper subword on of tt such that the values occurring in on form an interval (of 
course, oti might well be a singleton). 

Suppose that there are i < j such that on and aj overlap properly but are not 
equal. Then the elements of tt belonging to cither on or oej form a subword 
a whose values are the union of two overlapping intervals, hence an interval. 
By the maximality of either on or aj it must be the case that a — tt. Then 
on and aj with the elements common to on deleted form either a 12 or a 21 
decomposition of 7r. These cases are clearly mutually exclusive. 

Henceforth suppose that tt is neither 12-decomposable nor 21-dccomposable. 
Then the on form a partition of tt (i.e. any two are either equal or disjoint). 
The relative ordering of the on must be some absolutely irreducible permutation 
9, for otherwise we could pool some proper subset of the a's to form a coarser 
partition, contradicting the choice of each on as a maximal proper subword 
of tt whose values form an interval. Now reindex the distinct on and write 

TT = OL\OL<2 ■ ■ ■ Oik- 

Thus we have established the existence of a decomposition of the type claimed. 
To establish uniqueness, suppose that another decomposition of the same kind, 
say n — (3i(3 2 • • • Pm were given. We include here the assumption that the rela- 
tive ordering of Pi through [3 m forms an absolutely irreducible permutation. If 
Pi ^ ot\ then (3\ is a subword of a\ by the maximality of a\ . Now take the least j 
such that /3i/?2 ■ ■ ■ Pj contains a\. Then in fact we must have P1P2 ■ ■ ■ Pj = a\ for 
otherwise the values in a\ and pj form overlapping intervals, and so the values in 
P1P2 ■ ■ ■ Pj form an interval, contradicting the maximality of a\ (P1P2 ■ ■ ■ Pj 7^ w 
since n is neither 12 nor 21-decomposable). However, /3i/3 2 ■ ■ ■ Pj = ot\ contra- 
dicts the absolute irreducibility of the relative ordering of the /3's. So ot\ — Pi. 
But now the same argument implies that «2 = P2 and, inductively that in fact 
m = n and on — Pi for all i. m 

We now return to the analysis of F(321). Let tt e F(321) be given. Suppose 
that it is neither 12-decomposable nor 21-decomposable. By the proposition 
above, tt — ot\ot2 ■ ■ ■ for some subwords on whose values form intervals, and 
whose relative ordering forms an absolutely irreducible permutation. Since tt 
has some decomposition into subwords whose relative ordering avoids 321, and 
since the proof of the proposition above shows that the a's form the coarsest 
possible proper partition of tt into subwords whose values form intervals, it must 
be the case that the relative ordering of the «j avoids 321. Of course, we also 
have that each a% belongs to F(321). Conversely, given on in F(321), shifted 
to have relative order equal to some absolutely irreducible clement 9 of A(321), 
then, by the very definition of F(321), the permutation a\ct2 • ■ ■ ctk belongs to 
F(321). 

Let F(x) be the generating function for F(321), taken to have constant term 0, 
and Ai(t) be the univariate generating function for the absolutely indecompos- 
able members of A(321) of length greater than or equal to 3. From the above we 
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see that Ai(F(x)) is the generating function for the elements of F of length at 
least 3 which are neither 12-decomposable nor 21-decomposable. Let F + denote 
the generating function for the 12-indecomposable elements of F, and F_ that of 
the 21-indecomposable elements of F, again taken with constant term 0. Then 
F+F enumerates the 12-decomposable members of F while F_F enumerates 
the 21-decomposables. Further relations arise from the observation that a 12- 
indecomposable is either 21-decomposable or both 12- and 21-indecomposable, 
and similarly for minus indecomposables. We thereby obtain the system of 
equations: 



F = 
F+ = 
F_ = 



F+F- 

F^F 

F + F 



F F ~\~ A.j^'py* 

Ai rr \F) 
Ai rr \Fj . 



(F) 



Solving this system for F gives: 
F 2 + (A lrr (F) - 



1 + x)F + A irr (F) + x = 0. 



(7) 



We can use the work of the previous section to obtain a radical expansion of 
A ■ • 

1 — X — V — 3 X 2 — 2 X + 1 2 

A md {x) = ^Tl) x ■ 

Then substitution in (0) and elimination of radicals gives: 

F 6 + (-2 x + 3) F 4 + (-2 x - 1) F 3 + 
(~3x + 3 + x 2 )F 2 + (2x 2 - 1 - 2x) F + x + x 2 = 0. (8) 

The first few terms of the associated power series are: 

x + 2 x 2 + 6 a; 3 + 24 x 4 + 116 a; 5 + 625 x 6 + 3580 x 7 + 
21297 x s + 130084 x 9 + 810737 a; 10 + O (x 11 ) 

and the exponential constant governing the growth rate, is the reciprocal of the 
radius of convergence of this series. This radius is the least positive root of the 
discriminant of (||), which is an irreducible polynomial of degree 7. The value 
of the exponential constant is approximately 7.346751, compared to 4 for the 
underlying class A(321). 



5 Summary and Conclusions 

We began this paper with an explicitly constructed grammar to describe the 
skeletons of 321 avoiding permutations. In general there is a close connection 
between combinatorial classes with algebraic generating functions and unam- 
biguous context free languages. This connection can either be used, as here, 



15 



to provide an explicit enumeration of a class, or to provide a "soft" proof that 
the generating function of a class is in fact algebraic. The former approach 
has become much more attractive with the ready availability of symbolic alge- 
bra packages since the algebraic manipulations necessary to solve the equations 
arising from the grammar are undeniably tedious. The latter approach has been 
used in Q to provide algebraicity results for a family of pattern classes. It can 
also be used in the context of generating functions for generating trees, thereby 
generalizing a number of the theorems in |Q] about the existence of algebraic 
generating functions. It must be noted though that the results of that paper and 
similar results in Q provide much more explicit detail concerning the generating 
functions that they produce. 

One of the striking features of the equations for the various irreducible and 
indecomposable subsets of ^4(321) is their simplicity. In some sense then the 
cnumerative coincidences that we observed are not so startling, since there is a 
relatively limited supply of simple quadratic equations. Because of the simple 
form of these equations, one could apply Lagrange inversion to obtain explicit 
formulae for many of the coefficients, as is done for example in or indeed 
carry out detailed asymptotic analyses of these coefficients. 

We restricted ourselves to bivariate generating functions but the reader should 
note that the techniques employed can be naturally applied to produce other 
statistics of these permutations. For example, it would be a simple matter 
to produce, if one wished, a generating function A(x, y, z, w) where x marked 
total size, y marked left to right maxima, z marked the number of occurences 
of i(i + 1) among the left to right maxima, and w that number among the 
remaining elements. 

The class A(321) is the simplest pattern class, in terms of the patterns which 
it avoids, that contains infinitely many absolutely irreducibles. The techniques 
used in the preceding section to solve (in the sense of enumeration) the wreath 
fixed point equation: 

X = AlX 

apply, owing to proposition |J, completely generally to any base class A in which 
the absolutely irreducibles can be enumerated. In particular for a positive inte- 
ger n let D n be the class of permutations which "fractally have < n elements" . 
That is, they are comprised of at most n blocks, each of which is comprised of 
at most n blocks, each of which . . . Then D n is the solution of the fixed point 
equation: 

X = F n lX 

where F n is the class of permutations of size < n. Since F n is finite equation 
(0) is simply a polynomial and we obtain: 

Corollary 5 Each of the classes D n has an algebraic generating function. 
There is much further information to be gleaned from the representation of a 
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class as a subclass of D n when this is possible, and we hope to explore these 
matters in a future paper. 

One aspect of F(321) that has been notably omitted is a description in terms 
of minimal forbidden patterns. It appears that this set may be finite consisting 
of: 

42513, 35142, 41352, 362514, 531642 

but all that can be said with certainty at this point is that no further minimal 
forbidden patterns exist of length 12 or less. 
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